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Measurements are presented ofboundarv’ layers with embedded vortices and with 
film cooling for freestream velocities of 15, and 11 m/s. Measurements of a boundary 
layer with embedded vortex and and no film cooling, and of a boundary layer with film 
cooling but no vortex are presented for freestream velocity of 15 m s. Plots of total 
velocity, V, streamwise velocity, V.^, secondary (low vectors, total pressure, P^, and 
streamwise vorticity are presented for many of these test conditions. 

The results show that the embedded vortices completely dominate the (low field 
in boundary layers with film cooling. This is indicated from the plots of V, N ^ and P^ 
which show the effects of film cooling to be completely decimated in the vicinity of the 
vortex. 

In order to conduct this study, a five-hole pressure probe was calibrated for pitch 
and yaw. The probe was then used to measure five pressures associated with the flow. 
From these pressures, total velocity and the x, y. and z components of velocity were 
determined. 

A boundary layer profile was conducted to verify the calibration of the pressure 
probe, measurement procedures, and velocity computations. The results show expected 
boundaiy layer behavior with a small and component. 
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I. INTRODUCTION 



The increasing need for greater efTiciency in gas turbine engines has resulted in 
higher turbine inlet temperatures. Consequently, combustor liners and turbine blading 
are subjected to greater amounts of thermal stress, thermal fatigue, and creep. At 
present, gas turbines, such as those associated with militaiy' applications, have inlet 
temperatures as high as 1800 - 2000 degrees C (3270 - 3632 F) with pressures of 40 
atmospheres. 

Turbine parts may be protected from heat loads resulting from exposure to gas at 
high temperatures by using a coolant within turbine passages and along turbine 
surfaces. Convection cooling, impingement cooling, transpiration cooling, and film 
cooling are used for this purpose. Although it is possible for the cooling medium to be 
a substance such as liquid water or Freon - 12, most gas turbine arrangements utilize 
engine air bled off from the compressor and rerouted to the turbine nozzles and 
blading. With film cooling, compressor air is ejected from surfaces of blades and 
vanes. The film coolant then protects metal surfaces by forming a protective insulating 
film between the blades and the mainstream, and by acting as a heat sink. 

The flow through a turbine cascade is extremely complex. Efforts to analytically 
model the flow are successful only when the fluid is considered to be inviscid. When 
the effects of viscosity are included, the analysis is much less eflective [Ref 1]. Because 
of the difficulty involved in analytical representation of the viscous portions of cascade 
flows, much of the work in this area is experimental in nature. 

Flow visualization studies by Herzig, et al., [Ref 2] were among the first to show 
the complexity of flow within the cascade. Flow visualization studies by Langston, et 
al., [Ref 3] and Marchal and Sieverding, [Ref 4] also show the detailed development 
and complexity of the flow through the turbine cascade. More recently. Sieverding and 
Van Den Bosche, [Ref 5] have used color smoke-visualization to study the evolution of 
flow in cascades. 

Figure l.l from [Ref 6] shows the various secondary flows associated with flow 
in turbine cascades. As the inlet boundary layer approaches the blade's leading edge, a 
horseshoe vortex is formed. The point at which this formation occurs is the saddle 
point which is clearly shown in Figure 1.2 from [Ref 4]. One leg of the horseshoe 
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vortex moves initially near the side of the blade, and then through the blade passage 
where it is referred to as the passage vortex. In cascade flow, the passage vortex is 
composed of 1) fluid from the pressure side leg of the horseshoe vortex, 2) the 
crossflow from the endwall boundary layer, and 3) entrained fluid from the mainstream 
flow, [Ref 7]. As the passage vortex continues through the cascade, it is forced by the 
pressure gradient to the opposite side of the passage near the suction side of the 
adjacent blade. This shift is clearly shown in Figure 1.3 from [Ref 2]. The second leg 
of the horseshoe vortex follows the suction side of the blade and has a sense of 
rotation opposite to that of the pressure side vortex. This vortex moves away from the 
corner and is generally believed to be smaller in size than the passage vortex. As can 
be seen from Figure 1.4 from [Ref S] the passage vortex makes approximately one 
rotation as it passes through the cascade. In Figure 1.1 the number of revolutions has 
been exaggerated for clarity. 

The objective of this thesis is to study the effects of embedded vortices on a film 
cooled turbulent boundan,’ layer. In order to understand the effects that this complex 
flow field has upon heat transfer, vortex characteristics and their interaction with 
surrounding flow must be understood. In this study, a five-hole pressure probe is used 
to measure vortex characteristics. 

Extensive procedures for calibration and qualification of the five-hole pressure 
probe, and its use in measuring three-dimensional flows are first discussed. Results of 
six different test are then given: 1) a baseline measurement in a developing boundaiy 
layer. 2) a boundary’ layer with film cooling only at a freestream velocity of 15m s, 3) a 
boundary' layer with a single embedded vortex and without film cooling at a freestream 
velocity of 15 m s. 4) a boundary layer with embedded vortex and film cooling at a 
freestream velocity of 15 m. s, 5) a boundary layer with an embedded vortex and no 
film cooling at a freestream velocity of 20 m/s, and 6) a boundary layer with embedded 
vortex and film cooling at a freestream velocity of 11 m s. Cases 3, and 4 provide the 
most informative results of the study; cases 5 and 6 are incidental. 
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II. EXPERIMENTAL APPARATUS 



A. WIND TUNNEL 

The wind tunnel pictured in Figure 2.1 is now described. It is an open-circuit 
blower tunnel used to provide uniform flow at the nozzle exit. 

1. Description 

The wind tunnel is designated the NFS Shear Layer Research Facility (SLRF) 
and was built by Aerolab. It is designed to provide uniform flow with a minimum 
amount of turbulence intensity. It is designed with numerous pressure taps and four 
38 X 20.3 cm (15 x 8 in.) access ports along each of the of the side walls. The height of 
the top wall is adjustable to permit changes in the pressure gradient along the length of 
the test section. Additionally, the top wall contains numerous instrument ports for the 
measurement of various flow characteristics. 

The air speed through the tunnel can be adjusted from 5 m/s to 40 m, s. The 
blower exit slips into the inlet end of the wide-angle difiuser with 1.6 mm of clearance 
so that the fan is isolated from the body of the wind tunnel. The dilTuser inlet contains 
a filter and nozzle. The test section is 3.048 m (10 ft.) long and 0.6096 m (2 ft.) wide. 
The top is fabricated from Lexan sheet (4.76 mm thick), continuously sealed with 
neoprene along the edges. The tunnel's bottom wall consists of one 1.2192 m (4 ft.) 
long and three 0.6096 m (2 ft.) long removeable and replaceable sections. These 
sections are all 0.6096 m (2 ft.) wide and are sealed with "O" rings around the sides. 
Further discussions of the wind tunnel are contained in [Ref 9] and [Ref 10: p. 38]. 

2. Qualification and Performance 

Extensive qualification test of the Shear Layer Research Facility were 
conducted by Ligrani, [Ref 11]. Results show that the variation of total pressure at 
the exit plane of the nozzle is less than 0.4% at 26 m/s and 34 m s. Mean velocity 
varies less than 0.7% for the same mean freestream speeds. From five-hole pressure 
probe measurements, the velocity angle deviation is nowhere greater than about 0.6 
degrees at the nozzle exit plane. 

Profile measurements of the mean velocity and longitudinal turbulence 
intensity in the turbulent boundary layer developing at 20 m, s indicate normal, 
spanwise uniform behavior. For this qualification test, and all results which follow, the 



boundan' layer was tripped near the exit of the nozzle with a 1.5 mm high strip of 
tape. Total pressure measurements along the test section surface at the nozzle exit 
were uniform within 0.5% indicating spanwise uniform skin friction. 

Freestream turbulence intensity was measure to be 0.00085 (8.5 one - 
hundredths of one percent or .085 precent) at 20 m, s increasing to 0,00095 at 30 m, s. 

B. INJECTION SYSTEM 

Ordinarily, the injection system provides film coolant at temperatures above 
ambient. The freestream air is at ambient temperature; therefore, the heat transfer 
would then be in a direction opposite to that which occurs in gas turbines. The 
coolant is injected into the boundarx’ layer through a single row of injection holes. The 
injection holes are scaled relative to boundaiy layer thickness to be similar to those 
used in current turbine blade design. 

For the present tests, all injected air was at ambient temperatures. 

Injection system air is provided by an 71TD Ingersoll-Rand air compressor. The 
air is discharged from the compressor into three large storage tanks. As the schematic. 
Figure 2.2. shows the air flows from the storage tank through an adjustable regulator, 
a cut-olT valve, moisture separator, flow regulator, a Fisher and Portor rotometer (full 
scale 9.345E-3 m;s, 19.8 SCF.Vl, model 10.‘\3565A). The rotometer which controls the 
volumetric flow rate, discharges the film coolant through a diffuser and into the 
injection heat exchanger and plenum chamber. 

The heat exchanger and plenum chamber shown in Figure 2.3 is 0.305 x 0.508 x 
0.457 m (23 x 20 x 18 in.) and is constructed of 1.27 cm (1/2 in.) plexiglass. Injection 
air flows over three metal plates 0.381 x 0.508 m (15 x 20 in.). The two lower plates 
are covered by silicon rubber heaters, 0.381 x 0.483 m (15 x 19 in.), rated at 120 volts. 
The heaters are controlled through a type 136 Powerstat variable autotransformer. 
The top surface of the chamber contains 13 plexiglass injection tubes each being 8 cm 
(3.15 in.) in length and with an inner diameter of 0.95 cm (3 '8 in.). This corresponds 
to a length to diameter ratio of 8.42. The 13 injection holes discharge the coolant into 
the boundary' layer at a 30 degree angle. There is a three-diameter spanwise spacing 
between the center of each hole. 

Further discussion of the qualification and performance of the injection system 
can be found in [Ref 10: p.23]. 
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C. FIVE-HOLE PRESSURE PROBE 

Multi-hole pressure probes are invaluable in the investigation and measurement 
of complex, three-dimensional Hows. In particular, the five-hole probe is well suited 
for measurement of three mean velocity components in low speed incompressible flows. 

The five-hole pressure probe used to measure pressure in this study is 
manufactured by United Sensors and Control Corp. (drawing number 

D. '\-125-24-F-22-CD). The probe shown in Figure 2.4 is 0.6096 m (2 ft.) in overall 
length with a probe diameter of 0.318 cm (0.125 in.). It is constructed of corrosion- 
resistant, non-magnetic stainless steel. 

The five pressure holes are arranged in two different planes which intersect at the 
mutual hole, Pj. The probe tip is prismatic in geometiy, as shown in Figure 2.4 The 
centrally located Pj hole is normal to the freestream. The pitch plane consists of Pj, 
P^, and P^ pressure holes. While pressure holes Pj, P2, P3 constitute the yaw plane. 
The distance seperating P2 and P3 is 0.178 cm (0.070 in.), P^ and P3 are 0. 155 cm 
(0.061 in.) apart. The central hole Pj is 0.648 cm (0.255 in.) from the bottom of the 
probe tip. These distances were measured using a micrometer. 

For calibration, the probe was positioned in a manual traversing unit 
manufactured by United Sensor and Control Corp., Figure 2.5. The unit was modified 
by the user to include a compass rose with a radius of 15.24 cm (6 in.) and a range of 
yaw angles from -40 to -1-40 degrees (0.25 degree accuracy). Additionally, the manual 
traversing unit positions the probe vertically in the mainstream 0 - 30.48 cm (0 - 12 in.) 
range. The manual traversing unit is mounted on top of a spanwise, horizontal sled. 
The sled is 62.23 cm (24.5 in.) wide and is designed to set in place on top of the wind 
tunnel's side walls. The spanwise mounting sled is designed such that the probe can be 
positioned 10.16 cm (4 in.) either side of centerline in increments of 0.635 cm (0.25 in.). 
In addition to spanwise positioning, the spanwise sled could be rotated through a range 
of pitch angles from -15 to + 15 degrees with an accuracy of 0.5 degrees. 

After the completion of the probe calibration, an automated traversing 
mechanism. Figure 2.6, was used for probe positioning while measuring the pressures in 
the experimental test cases. The probe is fixed into the automated traversing 
mechanism in a position of zero yaw. The traversing mechanism has two degrees of 
movement which allows a thorough measurement of the flow field to be conducted. 
Both the spanwise and vertical traversing blocks are mounted on a 20-thread per inch 
drive screw and two ground steel, case-hardened steel guide; support shafts. Each drive 
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shaft is directly coupled to a SLO-SYN type Y1O92-FD310 stepping motor. The 
motors are controlled by a VIITAS Two-Axis Motion Controller, Figure 2.7 The 
stepping motors and the controller are manufactured by Superior Electric Company. 
The controller directs the movement the probe in both the spamvise amd vertical 
directions. The MITAS controller comes equipped with 2K bytes of memory and an 
.MC68000, 16-bit microprocessor which allows the user to program the start, stop, 
duration, speed, acceleration and deceleration of the stepping motors. 

D. DATA ACQUISITION SYSTEM 

The data acquisition system, shown in Figure 2.8 rapidly acquires the voltages 
associated with each pressure, converts each voltage to pressure. 

1. Transducers and Demodulators 

The probe is connected through reinforced plastic tubing to five Celesco model 
LCVR differential pressure transducers. These transducers have a designed pressure 
range of 0 -20 cm (0 - 7.85 in.) water and produce a 15 to 45 mV volt output signal. 
The transducer output signal is converted to a proportional DC signal by Celesco CD 
lOD carrier demodulators. Each demodulator has a maximum frequency response of 
-3dB at 500 Hz and a maximum out put noise of lOmV, peak to peak. Each 
transducer, carrier demodulator combination was calibrated against a Vleridian 1.27 cm 
(0.5 in.) horizontal manometer with an accuracy of 0.002 cm (.005 in.) of water, to give 
typical calibration of approximately l.O volt per inch of water differential pressure. 

2. Computers and Hardware 

A Hewlett-Packard 85 microcomputer was used to acquire and process data 
for the calibration of the pressure probe. Configured with 64K bytes of memorv' and a 
single magnetic tape cartridge drive, the HP-85 was used to collect, store, display, and 
print the majority of the data required during the course of probe calibration. 

For the measurement of the flow field, a Hewlett-Packard Series 300, Model 
9836S computer was dedicated to the data acquisition process. The HP 9836S is 
equipped with a MC68000, 8 MHz 16 32-bit processor. Dual 5-1,4 inch floppy disk 
drives, and IM bytes of memory. A HP 7470 two pen plotter was used for the graphic 
representation of data. 

Each transducer, 'carrier demodulator combination is connected directly to a 
HP 3498A extender which is controlled by a HP-3497A data acquisition'control unit. 
The HP-3497.A which provides precision measurement and process monitoring, is 
equipped with analog multiplexing and a digital voltmeter with l^V sensitivity. 



15 



Six software programs were developed for use during the thesis. PROCAL 
was developed for use with the HP-85 and was used during the probe calibration phase. 
PRSACQ, VEL, VELC, PLOT, VECTOR were developed for use with the HP 9S36S. 
PRSACQ is used to measure the pressures in the various flow fields. VEL and VELC 
are used to compute the velocity components. PLOT and VECTOR are used for the 
plotting of results. A thorough discussion of the programs requires an understanding 
of the calibration procedures and the velocity measurement techniques as discussed in 
Chapter Three. 
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III. CALIBRATION AND MEASUREMENT PROCEDURES 



Prior to using the five-hole pressure probe to measure flow velocities, it is 
necessarv' to calibrate the probe to determine the dependence of yaw. pitch, static, and 
total coefiicients on yaw and pitch angles. 

A. COORDINATE SYSTEM 

A right hand coordinate system, Figure 3.1 was established for use throughout 
the course of study. The X-axis is parallel to the streamwise direction and is positive in 
the downstream direction. The Y-axis is in the vertical plane and is positive from the 
wind tunnel's bottom wall. The Z-axis is in the vertical plane and is positive from the 
wind tunnel's bottom wall. The Z-axis is in the spanwise direction. The origin of the 
coordinate system is located on the centerline line of the bottom wall. 

Yaw, p. is defined as rotation about the Y-axis and was arbitrarily defined as 
positive when the direction of the flow resulted in pressure being greater than P 2 . 
This condition corresponds to positive component of velocity in Z direction. The pitch 
angle, a, is defined as rotation about the Z-axis and is defined as positive when P^ is 
greater than P^.This condition results in a positive component of velocity in the Y 
direction. 

B. EXPERIMENTAL PROCEDURES 

The five-hole pressure probe was calibrated using the method described by 
Treaster and Yocum. [Ref. 12). 

To perform the calibration, the probe is placed in the manual traversing unit 
which is mounted on top of the spanwise horizontal sled as described in Chapter Two. 
The probe is normal to the freestream when P 2 is equal to P 3 which gives p equal to 
zero. 

The probe was manually fi.xed at a predetermined yaw angle and then rotated 
through the pitch plane. At each point, the data acquisition system records the five 
pressures. Four pressure coefiicients are then calculated by the PROCAL program. 

Calibration of the probe was conducted over a range of yaw angles from -20 
degrees to +20 degrees in four-degree increments. The pitch angle was varied from -15 
degrees to +15 degrees in five-degree increments. This provided a cone of angles 
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which was suHicient for the Hows to be studied in this thesis. Calibration was 
performed in the NFS Shear Layer Tunnel at a freestream velocity of approximately 21 
m,'s (68.89 ft/sec). The probe was located 1.6 m (5.3ft) from the boundary layer trip. 
This equates to a Reynolds number of 1.97x10^ based on downstream distance. Static 
pressure was obtained from a static pressure tap on the tunnel side wall and total 
pressure from a Kiel probe inserted into the flow through the top wall. Static and total 
pressures are measured once per set of calibration points. Atmospheric pressure is 
used as the reference pressure. 

C. CALIBRATION COEFFICIENTS 

To use the probe for measurement of complex flows, it is necessary to determine 
the flow angles, a and P, and the local static and total pressures. This can be done by 
determining four non-dimensional pressure coeflicients over a range of angles in both 
the yaw and pitch planes. The four calibration coefficients are defined as: 



CPyaw= (P2-P3)(Pl"P) (eqn3.1) 

CPpitch = (P4-P5)'(Pl~P) (eqn3.2) 

CPtotal = (Pi ■ Ptotal)'(Pl ■ P) (®q^ 

^Pstatic ~ (P ~ Pstatic^'(Pl ~ P^ (eqn 3.4) 

where 

P = (P 2 + P 3 + P4+P5)/4 (eqn 3.5) 



To be of value in measuring a flow field, the calibration coefTicients must be a 
function of flow angle only, independent of velocity and repeatable. Figure 3.2 shows 
that the calibration is both independent of velocity and repeatable in the yaw plane. 
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As previously stated, probe calibration was conducted by fixing the yaw angle 
and varying the pitch angle. Theoretically, a calibration procedure of fixing the pitch 
angle and rotating through the yaw angles should provide identical results. However, 
in this study when the fixed pitch, vary yaw method was attempted, the results were not 
identical. The trend of Cpp-^^j^ vs. pitch angle showed significant scatter, and 
vs. pitch angle was not constant for each yaw angle. 

D. RESULTS 

Figure 3.3 shows the variation of with the yaw angles for various pitch 

angles. The response is nearly linear for all yaw angles. The results show that values 
of for various pitch angles collapse on top of each other for yaw angles of -12 

degrees to + 8 degrees. This means that is independent of pitch angle in this 

range. From -23 degrees to -15 degrees and from +8 degrees to +15 degrees, there 
are slight variations in the results indicating that has a slight dependency on the 

pitch angle. The fact that yaw is only slightly dependent on the pitch angle allows data 
to be more easily processed in determination of flow velocity. 

The variation of Cpp-^^j^ vs. pitch angle, Figures 3.4 and 3.5, show that this 
coefficient is dependent on both yaw and pitch angles. The trend of Cppj^^j^ vs. pitch 
angle is generally linear but there are variations from that linearity for each yaw angle, 
and unlike the yaw plane, these variations are not restricted to any particular region. 

Figures 3.4 and 3.5 show that the range of values for Cpp-^^j^ is small over the 
range of a shown compared to the variation of with p. [Ref. 12; pp. 27-28] 

attribtues this to the types of surfaces on which the holes in the pitch and yaw plane 
are connected. Large yaw angles result in one hole being nearly aligned with the flow. 
This hole senses a pressure nearly equal to the total pressure of the freestream. The 
other hole is then blocked from the freestream and, consequently, reads a pressure 
much less than the freestream static pressure. The holes in the pitch plane have a 
different response to variation in the pitch angle. When pitched, one hole reads a 
pressure which is near the total pressure of the freestream, but the second hole senses a 
pressure which is greater than the freestream static pressure. Thus, is a much 

larger number than Cpp^^j^. The small range of Cppjj^,j^ increases the scatter and 
uncertainty of pitch angle measurements. 

The plot of vs. pitch angle. Figure 3.6, shows that has weak 

dependency on yaw and pitch angles. 
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^Ptotal angle, Figure 3.7, indicates that for any given vaw angle, 

^Ptotal ‘s constant throughout the range of pitch angles. Here, also shows a 

weak dependency on yaw and pitch angles. 

E. INTERPOLATION AND APPLICATION 

After the probe has been calibrated and the operating characteristics of the probe 
are known, it is possible to determine the pitch angle, yaw angle, local static and total 
pressures for any flow field. 

The probe is positioned normal to the freestream, and at any location in the flow 
field, the five pressures can be measured. These five pressures are now used to 

calculate the experimental or local coefficients of yaw, pitch, static pressure, and total 
pressure. 

A fifth order polynomial was fitted to the average values of Cp The 

resulting polynomial computed using a FORTRAN program based on the leas^t-'squares 
method is; 

Papp = ■ ■ 1 58 - 7.36( M j ) + 0. 1 35( M j )2 + 0. 304( M j + 0.009{ M j )^ - 0.03 1 { M j )5 (eq n 3 .5) 

Here, P^pp is the approximate yaw angle, and Mj is the local coefficient of vaw. 

Knowing the approximate yaw angle and the local a computerized 

interpolation is performed to determine the pitch angle. Since the value of the pitch 
angle is dependent on the local Cppj^^^ and on the yaw angle, it is necessary^ to 
perform a double interpolation. Referring to Figure 3.8, the two yaw angles, Pj and 
p 2 ^ which bracket the approximate yaw angle are first determined. Next, the local 
^Ppitch, ^^2a ^*2b, located between known values of Cpp|j.^j^ from the 

calibration data for each yaw angle. In Figure 3.8, these values are designated Cjj 
^12, ^ 21 , and C 22 which correspond to pitch angles from the calibration data of ajj’ 

" 12 , «21, «22- bracketing pitch angle is determined by the following 

interpolation 

(Cii -M2g)/(Cii -C, 2 ) = (ajj -ajp);{ajj -aj^) (eqn 3.6) 

which rearranged gives 

“if - (ec,n3.7) 
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The other bracketing pitch angle « 2 f is found in a sinnlar manner which gives 

^ 2 f ~ ®21 ^^21 ^ 22^^^^21 ^^ 2 b ^ ^^21 ^ 22 ^^ J . 8 ) 

The pitch angle for the flow apis found through a second interpolation which results in 
the relationship 

«f= «if- («if- «2fX(Pi - Papp)XPi - P 2 )) (eqn 3.9) 

Because has a slight dependency on pitch, it is necessar>’ to refine the 

approximate yaw angle once the local pitch angle is known. This is done with a 
computerized interpolation routine which is ver\' similar to that used for the pitch 
angle. As shown in Figure 3.9, the pitch angles from the calibration data, ttj and 
which bracket the local pitch angle, apare first determined. The local is 

then located between the values of from the calibration data for ttj This results 

in C^j and C^j which correspond to yaw angles Pjj and Pj 2 - Then, by linear 
interpolation Pjpis found 

(C31 - Mjj,);(C3j - C2i)= (((Pi 1 - Pif);(Pi 1 - P12)) (eqn 3 . 10 ) 

or 

Plf=Pll-(Pll-Pl2)((C3i-Mn,)/(C3i-C2i)) (eqn 3.11) 

The same interpolation is done for U 2 which gives 

P2f~ P 2 I ~ ^P21 ~ P22X(^32 ~ ^^32 ~ ^22^^ 3.12) 

P->p and P 2 P are the yaw angles which bracket the yaw angle of the flow. The 
flows yaw angle, Pp can be found by a second interpolation which results in the 
relationship 



Pf p2f ^P2f PifX(®2 ®f)X®2 



(eqn 3.13) 
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Since local total and static pressures are not measured at each probe location, it 
is necessar\’ to perform another double interpolation to compute their values. From 
the calibration data, the values of for each of the bracketing yaw angles and 

pitch angles are known. In Figure 3.10, these points are designated Cjj, Cj 2 - C 2 j, 
and C-)-). An interpolation is performed to determine the Cp^Q^^j corresponding to the 
local pitch angle. These points are designated, Cj^and C 2 p The local M 3 , is 

found using the relationship 

(Cjf- M 3 );(C,f- = (pj - pf)/(pi - P 2 ) (eqn 3.14) 

which when rearranged as shown below gives the local Cp^Q^^j for the flow 

M 3 = Cjf- (Cif- C 2 f)((Pj - Pf) (P, - P 2 )) (eqn 3.15) 

The local M^ is found in a similar manner. Referring to Figure 3.11, the 

relationship for determining the local Cp is then given by 

.M^=C 3 f-(C 3 f-q^)((Pl-pp (Pj-P,)) (eqn 3.16) 

F, VELOCITY COMPONENT DETERMINATION 

Once the values of the local Cp^Q^^^j and have been determined, the local 

total and static pressure can be calculated and, subsequently, the total velocity at the 
probe tip can be determined. 

The defining relationships for Cpj.Q^gj and Cp^^-g^j^ can be rearranged to 
determine and 

Ptotal=Pi-(CPtotalXPl-P) (^^^-17) 

Pstatic” P~^^PstaticXPi “ P) 

By using Bernoulli's equation, the magnitude of the local total velocity is: 



Y -y 






(eqn 3.19) 



The three components of velocity can now be determined using the total velocity 
vector and the local pitch and yaw angles. Referring to Figure 3.1, these velocity 
components are given by, 

= VcosacosP (eqn 3.20) 



Vy = Vsina 



(eqn 3.21) 



= Vcosasinp 



(eqn 3.22) 



G. SPATIAL RESOLUTION CORRELATION 

The calculation of may be influenced by the local total velocity gradient. 
Corrections for this effect may be made using the following relationship 

Here is the uncorrected value ofV,, and V , is the value of V„ corrected for 
>0 > y* y 

spatial resolution. The value of ly used was slightly greater than the distance between 
P_j and P^ (0.155 cm): a value of 0.200 cm gave constant through the two- 
dimensional boundary layer. 

H. SOFTWARE 

Six programs were developed for use during this study. They are PROCAL, 
PRSACQ, VEL VELC, PLOT, and VECTOR. Each program is written in BASIC. 
PROCAL was used with the HP-85, all others were written for the HP-9S36S. All of 
the programs are listed in Appendix C. 

PROCAL is a BASIC program used for the calibration of the pressure probe. 
The program begins by computing the correction factor for random noise associated 



with each transducer. The user is then prompted to manually calibrate each transducer 
against a horizontal monometer. Static pressures are input utilizing a static pressure 
tap on the side wall and one of the transducers. Total pressure of the freestream is 
input using a Keil probe inserted into freestream and a transducer. The user is next 
prompted to input the ambient pressure in inches of mercur>’. After the user has 
positioned the probe at the desired angles of yaw and pitch, those angles are input into 
the program. 

The computer then acquires the five voltages from the data acquisition system, 
converts each voltage to a pressure in inches of water, and then calculates the yaw, 
pitch, total, and static coelTicients of pressure. Finally, the yaw angle, pitch angle, and 
the four coefficients are stored on a separate file and printed out utilzing the HP-85's 
internal printer. 

PRSACQ was used to acquire the pressures during the experiment. PRSACQ 
begins by prompting the user for the number spanwise and vertical data points and the 
resolution. A matrix of data point location is then computed. Next, the transducers 
are corrected for noise and calibrated against the manometer. Freestream static and 
total pressures are measured and ambient conditions are input. The program enters a 
loop which samples each pressure ten times per probe location. The local and 

Cppitch are computed. Probe position, and P| and P^gtal stored 

in a matrix. At the end of the data collection run. these values are read into a data file 
on a floppy disk. 

VEL is used to process the raw data acquired by PRSACQ. The data is first read 
into computer memoiy'. The program computes an approximate using a 

polynomial fit. Next, double interpolation subroutine is used to comput local 
A second interpolation subroutine is used to refine the value of Cp^.^^.. The values of 
the Cppj^^j^ and Cpy^^ are used in a third interpolation subroutine to compute the 
local and The total pressure is found from the definition of Cp^^jj^j. 

The velocity at the probe tip is computed, and the x, y, and z components of velocity 
are determined. Probe position (y and z coordinates), total pressure, total velocity, and 
V V and are stored in a matrix and then read into a data file. 

VELC corrects the component for spatial resolution. PLOT is used to 
generate graphs of streamwise velocity, total velocity, and total pressure. VECTOR is 
used to plot the secondary' flow vectors. 
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IV. EXPERIMENTAL RESULTS 



The study was conducted in three parts. The first was the measurement of 
baseline data consisting of boundary' layer profiles in a two-dimensional mean How 
field. The second was the measurement of the boundaiy layer characteristics at a 
freestream velocity of 15 m,s with 1) film cooling only. 2) with embedded vorte.x only, 
and 3) with both film cooling and embedded vortex. The third part was the 
measurement of boundarv' layer with vortex at 20 m/s, and measurement of boundarv' 
layer with vortex and film cooling at a freestream velocity of 11 m/s. 

A. BASELINE RESULTS 

The baseline boundary’ layer profiles were conducted at a freestream velocity of 
22 m/s. For these tests, the top wall of the wind tunnel was adjusted at 20 m, s 
freestream velocity so that a zero pressure gradient existed within 0.15 mm water along 
the length of the test section. Profiles were taken at three spanwise locations z — +2.54 
cm. z = 0 cm, and z = -2.54 cm. Figure 4.1 shows measurements of the streamwise 
velocity, V^. Figure 4.2 shows measurements of V^. Figures 4.3 and 4.4 show 
mesurements of the normal velocity, V^,. The results of the baseline measurements 
indicate behavior expected of a 2-D turbulent boundary’ layer, since the figures show 
mean profiles to be spanwise uniform for all three velocity components. Figures 4.2 
and 4.3 show that the and the corrected components are small and nearly zero 
as would be expected. Figure 4.4 is the plot of uncorrected for spatial resolution. 

B. I5M/S FREESTREAM VELOCITY RESULTS 

The investigation of the flow field was conducted by using the five-hole pressure 
probe to measure pressures at 800 points in a spanwise plane. Data was taken at 20 
different vertical locations, each having 40 spanwise locations. All 15 m s 
measurements were taken at a location of 1.49 meters from the boundary layer trip or 
0.39 meters from the heat transfer plate leading edge. The film cooling cases were 
conducted with injection air at 75% of full scale on the rotometer which corresponds 
to a blowing ratio (ratio of coolant to mass fluxes) of 0.50. 
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1. Boundar}’ layer with film cooling 

The results of boundary layer with film cooling only and no vortex are shown 
in Figures 4.5 to 4.8. Figure 4.5 shows contours of the streamwise velocity, Figure 4.6 
shows total velocity contours. Figure 4.7 is the plot for total pressure, and Figure 4.8 
shows the secondary flow vectors. Away from the wall, outside the boundaiy layer, the 
first three figures show spatially uniform behavior for V, and P^. Secondary flow 
vectors are ver>’ small eveiywhere in Figure 4.8. Near the wall deficits of V^. V, and 
Pq correspond to locations of the film cooling jets which are located about every 3.0 
cm from tunnel centerline. 

2. Boundary layer with vortex 

The vortex is generated by using a half-delta wing which is 3.0 cm high wdth 
7.5 cm chordand an angle of 18°. It is identical to vortex generator # 2 described by 
Joseph, [Ref. 10: p.76]. The vortex generator w'as located at a spanwise location of 
z = 4.79 cm (note that the direction of +z in [Ref. 10: p.75] is reversed). Figure 4.9 
shows the streamwise velocity results. Figure 4.10 is a plot of the total velocity 
contours. Figure 4.11 is a contour plot of total pressure, and Figure 4.12 shows the 
secondary flow vectors. The contour plots for V^,V, and P^ show significant deficits 
caused by the generator wake which is rolled up with the vortex. The center of the 
wake is located at y=3.3 cm and z=-3.05 cm. Figure 4.12 show's that the vortex 
center is located near the same location, with overall characteristics similar to a 
Rankine vortex. Figure 4.13 shows the streamwise vorticity contours, w'here the 
vorticity is calculated using 

= {dY^ldy)- {dY^ldz) (eqn 4.1) 

The vorticity is largest near the vortex center as expected. The circulation of the 
vortex is estimated to be 0.2708 m^s using the equation 

r = JWj^dA. (eqn 4.2) 

3. Boundary' layer with vortex and film cooling 

Figure 4.14 show's the streamwise velocity contours for this case. Figure 4.15 is 
the total velocity contour plots. Figure 4.16 shows the results for total pressure. 
Figure 4.17 show's the secondary flow' vectors, and Figure 4.18 is a contour plot of 
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vorticity. Figure 4.15, 4-16, and 4-17 show that the deficits for V.^, V, and from the 
film cooling are no longer present near the vorte.x. This result shows that the effects of 
film cooling are decimated by the vortex. The effect of the vortex is particularly 
pronounced near its downwash side. As for the previous case, Figure 4.18 shows that 
vorticity is again highest near the vortex center. Figure 4.19 from [Ref 10: p.99] shows 
that the results of this case are consistent with previous work by Joseph which shows 
high heat transfer rates on the downwash side of the vortex and low heat transfer rates 
on the upwash side of the vortex. The circulation for this case is estimated to be 
0.2708 m^ s. 

C. 20 M/S AND 10 M/S FREESTREAM VELOCITY RESULTS 

Figure 4.20 is a plot of streamwise velocity for the case of an embedded vortex 
with no film cooling with a freestream velocity of 20 m,'s. Figure 4.21 is a plot of the 
secondarx’ flow vectors for the same case. Figure 4.22 and Figure 4.23 are for a 
boundary layer with embedded vortex and film cooling at 1 1 m/s freestream velocity. 
The blowing ratio for the 20 m's case is 0.38, and the blowing ratio at 1 1 m, s was 0.68. 
The results for both cases show trends which are similar to those discussed above. The 
results at a freestream velocity of 11 m s are believed to be less reliable because of 
disturbances at the inlet of the tunnel which propagated to the test section during the 
time period the data was acquired. 
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V. SUMMARY AND CONCLUSIONS 



Measurements are presented of boundary layers with embedded vortices and with 
film cooling for freestream velocities of 15, and 11 m/s. Measurements of a boundary 
layer with embedded vortex and and no film cooling, and of a boundary layer with film 
cooling but no vortex are presented for freestream velocity of 15 m s. Plots of total 
velocity, V, streamwise velocity, V,^, secondary flow vectors, total pressure. P^, and 
streamwise vorticity are presented for many of these test conditions. 

The results show that the embedded vortices completely dominate the flow field 
in boundary layers with film cooling. This is indicated from the plots of V, and P^ 
which show the eflects of film cooling to be completely decimated in the vicinity of the 
vortex. This result is consistent with the heat tranfer results of Joseph, [Ref 10: p.54], 
which shows a localized hot spot at the wall near the same location. Future film 
cooling injection-hole arrangements in turbine blades must be designed to compensate 
for hot spots due to the vortices. 

In order to conduct this study, a five-hole pressure probe was calibrated for pitch 
and yaw. The probe was then used to measure five pressures associated with the flow. 
From these pressures, total velocity and the x, y, and z components of velocity were 
determined. 

A boundaiy layer profile was conducted to verify the calibration of the pressure 
probe, measurement procedures, and velocity computations. The results show expected 
boundaiy layer behavior with a small and component. 

It is recommended that flow visulization study of the interaction between the 
vortex and film cooling be conducted to enhance the understanding of this complex 
phenomena. 
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APPENDIX A 
FIGURES 




Figure 1.1 Endwall secondary’ flows. 
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y 2 =160 m/ 5 



1. Saddle point 

2. Separation line of inlet end wall boundary layer 
3a. End wall separation line of suction side lee of 

lead mu ed<ic vortex 

3b. End wall separation line of ‘'pressure side le'j 
of leading edee vortex" and or passage vortex 

4 . Start of suction side end wall corner vortex 

5. Low enerizy region 



Figure 1.2 Smoke and oil visulization for rotor blade. 
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Figure 1.3 Secondary flow deflection. 
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Figure 1. 




Rotation of passage and corner vortices in turbine cascade. 



32 



I 




Figure 2.1 P! iOU'^^raphs of wind tunnel. 




Figure 2.2 Schematic of injection air flow. 
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Figure 2.3 Schematic of injection plenum. 
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Figure 2.4 Photographs of pressure probe. 
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Figure 2.5 Manual traversing device. 
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Figure 2.6 Automated traversing mechanism. 
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Figure 2.7 Photographs of Two-Axis Motion Controller. 
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Figure 2.8 Photographs of' data acquistion system. 
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Figure 3.1 Coordinate system for flow measurement. 




Figure 3.2 Repeatabilitv of calibration results at dilTerent 
frcestream velocities for vs. yaw angles. 
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CP YAW VS. YAW 




Figure 3.3 Probe calibration. vs. yaw angles, 

freestream velocity or-?! m/s. 
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Figure 3.4 Probe calibration. Cp ^ vs. pitch anelc for -20 » to -4<> vaw, 

freestream v6iydity of 21 m, s“ 
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PITCH ANGLE 




Figure 3.5 Probe calibration. Cp^:,„u vs. pitch anelefor O® to + 20° vaw, 

freestrearn v6icrdity of 21 rn/s“ 
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CP STATIC VS. PITCH 




Figure 3.6 Probe calibration, vs. pitch angles, freestream velocity of 21 m/s. 
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CP TOTAL VS. PITCH 
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Figure 3.7 Probe calibration, vs. pitch angles, freestream velocity of 21 m/s 
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Figure 3.8 Interpolation for pitch angles. 




Figure 3.9 Interpolation for yaw angles. 
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Figure 3.10 Interpolation for local CpjQ^gj. 
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Figure 4.1 Baseline boundary layer results for streamwise velocity, V 




Figure 4.2 Baseline boundary" layer results for spanwise velocity, V . 
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Figure 4.3 



Baseline boundary’ 



layer results for corrected V^,. 
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Figure 4.4 



Baseline boundary layer results for uncorrected V 
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Figure 4.5 Streamwise velocity for boundarv layer without embedded vortex, 
75% film coolmg, freestream velocity 15 m/s. 
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Figure 4.6 Total velocity for boundary layer without embedded vortex 
75% film co'oling, freestream velocity 15 m,'s. 
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Figure 4.7 Total pressure for boundarv’ layer without embedded vortex 
75% film cooling, frcestream velocity 15 m,'s. 
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Figure 4.8 Secondary flow vectors for boundaiy layer without embedded vortex, 
75% film cooling, freedtream velocity 15 m/s. 
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Figure 4.9 Streamwise velocity for boundary layer with embedded vortex, 
without film coolmg, freestream velocity 15 m;s. 
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Figure 4.10 Total velocity for boundar>^ layer with embedded vortex, without 
film cooling, freestream velocity 15 m/s. 
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Fieure 4.11 Total pressure for boundary layer with embedded vortex, 
without film cooling, freestream velocity lo m/s. 
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Figure 4.12 Secondary' flow vectors for boundary’ laver with embedded vortex, 
without film cooling, (rcestream velocity 15 m/s. 
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Figure 4.13 Vorticity contours for boundar>’ layer with embedded vortex, 
without film cooling, freestrearh velocity 15 m/s. 
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Figure 4.14 Streamvvise velocity for boundarv layer with embedded vortex 
and film cooling," freestream velocity 15 m/s. 
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Figure 4.15 Total velocity for boundary' laver with embedded vortex 
and film cooling, freestrearh velocity 15 m/s. 
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Figure 4.16 Total pressure for boundary layer with embedded vortex 
and film cooling, frcestream velocity 15 m/s. 
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Figure 4.17 Secondar>’ flow vectors for boundary layer with embedded vortex 
and film cooling, (rcestream velocity 15 m/s. 
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Figure 4.18 Vorticity contours for boundary’ laver with embedded vortex 
and film cooling, freestream velocity 15 m/s. 
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Figure 4.19 Stanton number ratio for vortex at z= -4.79cm with film cooling. 
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Figure 4.20 Streamwise velocity for embedded vortex 
without film cooling , freestream velocity 20 m/s. 
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SECONDARY FLON VECTORS 
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Figure 4.21 Secondar>' flow vectors for embedded vortex 
without filmcooling, freestream velocity 20 m;s. 
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Figure 4.22 Streanuvise velocity for embedded vortex 
with film cooling. iVecstream velocity 1 1 m/s. 
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Figure 4.23 Secondarv flow vectors for embedded vortex 
with film cooling, freestream velocity 11 m/s. 
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APPENDIX C 
SOFTWARE 



4 ;■% 


REM 


P.ROGKAM PROCAL 


2 0 


REM 


THIS PROGRAM COMPUTES 


30 




THE COEFFICIENTS OF VAN 


40 


REM 


AND PITCH FOR THE FIUE 


5 0 


REM 


HOLE PRESSURE PROBE 


60 


REM 




701 


ReM 


DAVID EVANS NOV 1586 


801 


REM 




901 


REM 


VARIABLE NAMES 


i 00 


REM 




• 1 ;7« 


REM 


E<I> IS THE VOLTAGE 


1 w 01 


REM 


READ FROM THE DATA 


138 


REM 


ACOUISTION SYSTEM. 


140 


REM 




150 


REM 


A < I > I S THE C 0 N V E R S 1 0 N 


160 


REM 


FROM VOLTAGE TO INCHES 




155 


REM OF WATER . 


1 70 


REM 




1 SO 


REM 


P-:;r> WHERE 1 = 1 TO 5 


1 901 


REM 


ARE THE PRESSURES 


2 01 0 


REM 


FROM THE 5 HOLES 


2 1 0 


REM 


OF THE PRESSURE PROBE 


2201 


REM 




230 


REM 


P6=STATIC PRESSURE 


240 


REM 


P7=T0TAL PRESSURE 


250 


REM 


P V 6> = AVERAGE PRESSURE 


260 


REM 


Dl= P2-P3 


270 


REM 


D2=P4-P5 


2S0 


REM 


D3=P1-P';6> 


290 


REM 


Cl= YAW COEFFICIENT 


30101 


REM 


C2 = PITCH COEFFICIENT 


310 


REM 


C3=C0EFFICIEHT OF TOTAL 


3201 


REM 


PRESSURE 


330 


REM 


C4=C0EFFICIENT OF STATIC 
PRESSURE 


348 


REM 




350 


DIM 


C 1 < 1 00 > C2 <: 1 00 > .. C3 ( 1 00 ) 


360 


DIM 


C4 < 1 00 > w-<: 1 < 1 00 > .. Y <100 ;• 


370 


REM 




380 


DI3P 


•'ENTER DATE> <MMDDYY> " 


390 


IHPU 


T Ml 


480 


DISP 


"ENTER TIME,<HHMM>" 


410 


INPU 


T M2 


420 


REM 




430 


PRINT Ml.. .. M2 


440 


REM 





45i3 CREATE ‘'CaLDfll".2 
460 RbSIGH4 1 TO "CALDfil" 
470 ! PRINT# 1 Ml.. M2 
430 REM 
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4y0 REM iHIS SECTION COMPUTES 
500 REM THE CONNER'S ION FhCTOR 
510 RtM FOR VOLTAGE TO INCHES 
520 REM OF NRTER.THE USER MUST 
530 REM INPUT THE MANOMETER 
540 REM 

550 REM THIS SECTION COMPUTES 

tl5 REM THE ZERO OFFSET' COR 
bib REM RECTI ON FOR THE 
T R A N S ~ 

560 REM THE ZERO OFFSET COR- 

622 REM RECTI ON FOR THE TRA 
NS- 623 REM 
570 REM 

5S0 DISP “COMPUTE ZERO OFFSET “ 

5 -3 0 D I S P " C 0 R R E C T I 0 N . “ 

600 DISP " ! ! ! DISCONNECT ALL 

TUBING TO THE TRA 
NS DUCERSi!!” 

610 DISP " HIT CONT. WHEN READY" 
b 20 PA U 3 E 
630 CLEAR 

640 FOR 1=200 TO 204 

6 5 0! Z 2 = 0 

b 6 0 F U R • J — 1 T U 5 0 

670 OUTPUT 703 ; " A I " I " UT 1 “ 

580 ENTER 703 ; X2 
630 Z 2 = Z 2 "S- X 2 
700 NEXT J 
710 Z3=Z2.''50 
720 Z4<' 1-133 >=Z3 
730 NEXT I 
740 REM 
750 REM 

760 REM ;t:lc;<c:4:;t':4:t:LOOP#l:*:*:4:*:4;;*:4.1::|: 
770 REM 

7S0 FOR 1=200 TO 204 
730 11=1-133 

800 DISP “TRANSDUCER CALIBRATION 
ENTER THE MANOMETE 
R PRESSURE OH 

H20> ■' 

S10 DISP "ENTER THE PRESSURE" 

820 DISP "TRANSDDUCER NR . " , 1 1 
830 INPUT H2 
340 Z0=0 
S50 REM 

860 REM YYYY:«cL00P#2lc:*::*::t:*4:*:C:^;f 

870 FOR J=1 TO 50 

880 OUTPUT 709 ; " A I" ; I ; " VT 1 " 

330 ENTER 703 X 
300 Z0=Z0+X 
310 NEXT J 
320 REM 
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S* 3 0 E 1 = Z ti • 5 0 

340 FK I-139>=H2^' <E1-Z4‘:; I-133> > 
350 PR I NT " R " I .. " = " R ( I - i 35 :■ 

PRm cusp "IF H<;i> IS IJHSfiT 
PRESS 8 . " 

370 INPUT N2 

380 IF N2 = 8 THEN GOTCi 300 
330 NEhT I 
1000 REM 

1 0 1 0 R E M 1 1 4 ; -H 4c 4c 1 1 % :4c ;4c t :4: :4: 't 
1020 REM 

1030 REM ENTER THE VRLUES OF 
1040 REM STRTIC AND TOTAL 
1050 REM PRESSURE. 

1 0 6 0 Rl E M 

1070 CISP "INPUT STATIC PRESSURE 

.'TWrurc: nr uon •" 



1 080 
1 030 
1 100 
1110 
1 120 
1 130 
1 1 40 
1 1 0 
1 160 
1 170 
1 i 30 
1 1 30 
1200 
1210 
1220 
1230 
1240 
1250 



1260 
1270 
1230 
1230 
1300 
1310 
1320 
1330 
1340 
1350 
1 360 
1370 



DISP "USE TRANSDUCER NR . 1 " 
DISP "AND CHANNEL HR. 2G0" 
REM 

DISP "HIT CONTINUE WHEN" 
DISP "P STATIC PROPERLY" 
DISP "HOOKED UP." 

PAUSE 
L E H R 
B6=0 

FOR J=1 TO 50 
REN 

OUTPUT 70S .: " A I " .; 208 ; " '...'T 1 " 

ENTER 703 .: X6 

B6=36+X6 

NEXT J 

P6=B6''50 

REM 

DISP "INPUT THE TOTAL PRES- 
SURE.'; INCHES H2 

0 > " 

DISP " AND CHANNEL NR. 200“ 
DISP " " 

DISP " HIT CONTINUE WHEN" 
DISP " P TOTAL PROPERLY" 
DISP " HOOKED UP." 

PAUSE 
CLEAR 
REM 
B7 = 0 

FOR J=1 TO 50 

OUTPUT 703 .; " A I " .; 2O0.: " VT 1 " 
ENTER 703 .■ X7 
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1380 97=B7+K7 
1390 WEKT J 
1400 P 7 — E! 7 5 0 
1410 REM 
1420 REM 

1430 REM ENTER H M E; I E N T C G N D I T I o H 
1440 REM 

1450 DISP "ENTER PRMB < I N . HG > " 
1460 INPUT fil 
1470 REM 
1480 T=0 

1490 F 0 R J = 1 TO 5 0 
1 500 OUTPUT 709 ; " ft I" 1 09 " VT I " 
1510 ENTER 709 T1 
1520 T=T+T! 

1530 NEKT J 
1540 T2 = Tx'50 
1550 GOSUB 3010 
1560 REM 

1570 REN DEL P COMPUTED FROM 

1580 REM PTOTftL-PSTftTIC 

1590 REM 

1600 C=P7— P6 

1610 REM 

1620 REM 

1630 REM CONVERSION TO 31 UNITS 
1 640 ft 1 =ft 1 ;t:3385 . 82 



1650 C=Ct248.7 

1669 T4=T3+273.15 

1670 Rl=ftl/ <287;*:T4:- 

1630 Ul = <2:t:C/Rl 5 

1690 REM 

1700 REM 

1710 REM ENTER THE BIG LOOP 
1720 REM FOR READING PRESSURES 
1730 REM AND COMPUTING COEFF . 
1740 REM 
1750 REM 

1760 FOR K=1 TO 100 
1770 K1=K 

1730 REM ENTER THE VALUE OF THE 
1790 REM PITCH ANGLE AND YAW. 
1300 REM 

1310 DISP " ENTER THE VALUE OF 

THE PITCH ANGLE" 

1320 INPUT XKK) 

1330 REM 

1340 DISP "ENTER THE VALUE OF" 
1350 DISP "THE YAW ANGLE." 
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1360 INPUT V-aO 

1370 DI3P '• PITCH RNGLE ENTERED 

I S " . X 1 C K > 

1330 DISP "Yaw PNGLE ENTERED IS 

" .. Y < K .) 

1390 DISP " 

1900 DISP " DO YOU WhNT TO CHRNG 
E PITCH OR YRW IF 

SO PRE SS 1 " 

1910 INPUT N4 

1920 IF N4=l THEN 1800 

1930 REM THIS LOOP ACQUIRES EACH 

1940 REM PRESSURE 50 TIMES AND 

1950 REM TAKES THE AOERAGE . 

I960 REM 

1970 DISP "HIT 'CONTINUE' FOR" 
1930 DISP "DATA ACQUISITION." 
1990 PAUSE 
2000 CLEAR 
2010 REM 

2020 REM :t:;t-T:J::t:tL00P#3.tTlc:|c:»:**;«c;t; 
2030 FOR 1=200 TO 204 
2040 Z1=0 

2050 REM TTt:«c:*;T:LOOP#4T:4::<c:|::t;;|::t:T:t: 

L 0 6 0 FOR U = 1 i U 50 

2070 OUTPUT 709 " A I " ; I " OT 1 " 

2030 ENTER 709 X 

2090 Z1=Z1+X 

2100 NEXT J 

2110 REM 

2120 X=Z1P50 

2130 E<I-199>=X 

2 140 P < I - 1 99 > =A I - 1 99 > TE I - 1 99 ) - 
Z4< 1-199 :■ 

2150 NEXT I 

2 1 6 0 R E M % t ■*: 1 1 1 1 1 1 1 « 1 1 1 1 :<c .t jf. 1 4 1 1 
2170 REM 
2180 P<6>=0 



2 1 90 


FOR 1=2 TO 5 


2200 


P (5 > =p < I > + P < 


2210 


NEXT I 


2220 


P<:6>=P<6>/4 


2230 


D1=P<2;--P<3) 


2240 


D2 = P':4>-P(5) 


2250 


D3=P<; 1 >-P<6> 


2260 


Cl <K>=D1/D3 


2270 


C2(K.':'=D2PD3 


22S0 


C 3 '1 K > = < P < 11'- 


2290 


C 4 K > = 'i P ( 6 ) - 
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2300 


R£M 




2310 


PR I NT# 1 ; K 1 C K ^ .• V < K ) .. C 1 < K > 




. C2<K) 


' . C3 < K > .. C4 f > 


2320 


REM 




2330 


REM 




2340 


PRINT 


'•PITCH ANGLE IS%XKK 


2350 


PRINT 


"VfiW ANGLE IS".. Y<K;* 


2360 


PRINT 




2370 


PRINT 


" P 1 = " , P 1 > 


23S0 


PRINT 


" P 2 = " .. P < 2 > 


2390 


PRINT 


"P3=" '.. P<3> 


2400 


PRINT 


"P4=" .. P<4> 


2410 


PRINT 


"P5='‘ .. P<5> 


2420 


PRINT 




2430 


PRINT 




2440 


DISP •' 


LAST POINT??" 


2450 


OISP " 


IF SO PRESS 1" 


2460 


INPUT 


N1 


2470 


IF Nl = 


:1 THEN bOTO 2500 


2480 


NEXT K 




2490 


REM t t M t :4c t :4c 1 1 1 :4c :4c :4c :4c :4c :4c :4c 1 1 1 




1705 




2500 


REM 




2510 


RSSIGH# 1 TO :4: 


2520 


REM 




2530 


PRINT 




2540 


PRINT 


" :4: :4: :4: :4c :4c :4: :4c :4c 4c T ^c ;4: :4c :4c :4c :4c :4: :4c 




:4:4;4:#::4:4.:4:;4:" 


2550 


PRINT 




2560 


PRINT 




2570 


PRINT 


"FIVE HOLE PRESSURE 






PROBE CAL I BRAT 




ION. " 




2580 


PRINT 




2590 


PRINT 




2600 


PRINT 


"DATE OF RUN IS".. Ml 


2610 


PRINT 


"TIME OF RUN IS".M2 


2620 


PRINT 




2630 


PRINT 




2640 


PRINT 


"DENSITY-^KG/MS^ " 


2650 


PRINT 


R1 


2660 


PRINT 


"VELOCITY<M/S>" 


2670 


PRINT 


U1 


2680 


PRINT 


"PAMBCN/ME) " 


2690 


PRINT 


Hi 


2700 


PRINT 


"TAMB<C> " 


2710 


PRINT 


T3 


2720 


PRINT 


"P STATIC <IH. H20;-" 


2730 


PRINT 


P6 


2740 


PRINT 


"P TOTAL = (. IN. H20 J" 


2750 


PRINT 


P7 


27 60 


PRINT 




2770 


PRINT 


"PIT TAW CPI C 




P2" 
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2 7 2 0 R E ri 1 1 1 1 L 0 0 P # 5 >: 1 1 1 1 1 #. t 
27-30 POR 1 = 1 TO K1 
2S00 P-R I NT US I NG 23 1 O ; X 1 < I > .. Y •' 
J > .• C i i I ) t C2 ( I > 

2310 I MhGE mod .. 2:': .. mod .. 2X .. MD . DDDD 
.2 X.. MD.DDDD 

2S20 NEXT I 

2 3 3 0 R E 1 1 t 1 1 1 1 1 :4c 1 1 % :♦: % % :4c :ic % 

2840 PRINT 
2358 PRINT 

2860 PRINT "PIT YflW CP3 CP 
4 

2 3 7 0 R E r-1 ;4. % 1 1 :4: $ L 0 0 P # 6 1 1 :¥ 4: t X % .f. :t: 
2380 FOE I=i TO K1 
2338 PRINT USING 2900 X1<I>,Y^ 

I > .. C3 < I ;> .. C4 < I ) 

2300 I MfiGE MOD , 2X . MOD , 2X , MD . DDDD 
,2 X.. MD.DDDD 

2310 NEXT I 

2320 REM :*::t:4::4c4::4::4::^t4:.l::^;4:4:*:4:.4:.4:*.ic:4:4;4c 
t. 

2370 PRINT 

2340 DISP " THRNKS FOR USING ME" 
2350 DISP " HOPE YOU FOUND WHfiT" 
2960 DISP "YOU WERE LOOKING FDR" 
2370 DISP " REGRRDS" 

2330 DISP " HP-85B" 

2330 END 
7000 REM 

3010 REM MOLTfiGE TO TEMPERATURE 
3020 REM CONVERSION 
3030 REM 

3040 REM •..'<MV> TO TcC> 

3050 REM 

3060 El=T2:*cl000 

3070 T3 = 26 . 573:4:E1-1 . 936879:i:E 1 4;E 1 
+ .99785t.El*El*El- ,261277.t.El 
:4:Ei:4:EltEi 
3030 RETURN 
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20 

30 

43 

50 

EO 

73 

30 

GO 

130 

i LO 

i:o 

1-30 

140 

150 

ISO 

170 

iCQ 

ISO 

200 

710 

230 

230 

240 

250 

2B0 

270 

2S'i‘ 

290 

300 

310 

320 

330 

340 

3G0 

3S0 

370 

380 

330 

400 

410 

420 

430 

440 

450 

4B0 

47Q 

430 

430 

530 

■319 

i,; 2 0 



Rt.-l PROGRAM rRSACQ 

REM THI'3 PROGRAM ACQUIRES THE PRESSURES FROM 

REM A FIvE HOLE FRE2SURE PROEE AMO COMPUTES CRVAU .amO 

REM OBTAIilED FROM THE PROGRAM "FROCAL" 

REM 

REM AUTHOR DACE EUANS , JAMUAR/ , FESRUARr , 

REM 

REM VARIABLE rJAMES 
REiM 

REM Ev I ) 12 THE UOLTAGE READ FROM THE DATA ACQUiSTIOM S'tSTEM 
REM H(I) IS THE COMUEHSIOM FPCM UOLTnC-E TO IMChES QF WATER 
REr5 P(I) IS THE PRESSURE READ FROM Ti^E PIUE HOLE PRESSURE PROBE 
REM 1 = 1 TO 5 AND IN EACH OF THE ABOUE: IS ASSOCIATED WITH THE 
REM PRESSURES FROM THE FIVE HOLE PROBE 
REM 

REM FS=STATIC PRESSURE QF THE FREE3TREAM 

REM P7=T0TAL pressure OF THE FPEESTREAM 

PEM P<S;=AUERAGE PRESSURE 

REM Di=P2-r3 

REM D2=F4-PS 

PEM D3=Pl-F(b) 

REM 

REM C1=GAlI 8RATED YAW COEEFICIENT 
REM G2=CALI9RATE0 PIT’CH COEFFICIENT 
REM C3=CnLIBRATED TOTAL COEFFICIENT 
REM C4-CALlf?RATED STATIC COEFFICIENT 
Run 

REM Ml^MEA.SURED YAW COEFFICIENT 
REM M2 =MEaSURED PITCH COEFFICIENT 
REM M5--=MEASURED TOTAL COEFFICIENT 
PEM rt4=MEASUH£D STATIC COEFFICIENT 
PEM 

REM V=THG LuCAl TOTAL VELOCITY VECTOR 
REM V1 = THE >: COMPONENT OF VELOCITY 
REM U2=THE Y COMPONENT OF VELOCITY 
REM VZ=THE 2 COMPONENT OF VELOCITY 
REM 

REM X1=THE PITCH ANGLE FROM THE CAUieRATIuN DATA 
REM Y=THE YAW ANGLE FROM THE CALIBRATION OATn 
REM 

REM X4=THE PITCH ANGLE OF THE FLOW 
REM 

REM V1=THE VERTICAL LuCATrC'N OF THE PROEE TIP 
REM Z0=THE 5PANWISE LOCATION OF THE PROBE TIP 
REM Y3=THE YAW ANGLE OF THE FLOW 
REM 

REM A1=AM8IENT PRESSURE 

PEM T4=ABS0LUTE ambient TEMP'; FELVIN ) 

REM Rl=DENSITY (KG/rr3) 

sTM UB=VElDCITY of THE FREESTC'EAM 

REM 
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530 

540 

S50 

5B3 

570 

580 

530 

£00 

6*0 

520 

530 

840 

650 

BBO 

B73 

6S0 

S33 

700 

7iO 

'z - T/s 

759 

743 

750 

750 

770 

780 

730 

S00 

813 

373 

833 

840 

SSG 

8G3 

379 

889 

see- 

goo 

513 

920 

933 

340 

950 

360 

361 

370 

983 

393 

1C0( 

101 



0 1 n Y ( 1 390 ) , Z < 1 000 ' . -M i \ 1 098 • , m2 ( 1 0O-3 ; 

DIM Pl( lOOG),PBi 13005 

«EM 

REM 

IMP-JI ■•EMTER DnTF. . ( MMDDV V ;• “ , K2 

INPDT -ENTER I In£' . ; HNMM ; % n3 

PR INIER la 701 

PRINT -O^iTE 01- RLiN IS“,N8 

PR 1 NT r I ME OF p'JH I S , NS 

REM 

ir^ror -ENTER PIS 5F,ir;wr3E“ . M3 

INPUT "ErnER PTS VERTICAL" .N3 

REM N3 MUST BE AN EUEM iriTEGEP 

INPUT "ENTER SPnNUilSE .RESOLU nCN( iN , Z4 

INPUT “ENTER vertical RE50LUTI0N( IN > " , y4 

INPUT •■INITIAL 2 UN;", 2 3 

INPUT "INITIAL VvIN ",Y3 

PRINT "PTS SPANwISE" ,n3 

PRINT PTS VERTICAL" ,N3 

PRINT 3PANW I SE PE30LUT I ON : 1 M ; " . Z 4 

PRINT ••'VERTICAL RESCLUTION <IN)",Y4 

•PFINT "INITIAL ZUN.'".Z3 

PRINT "INITIAL yUM)“,Y3 

REM 

z: 1 >=Z3 
Y-' 1 ) = -Y3 
M3=r;3/Z 

FOR 17=1 TO N9 
TG=I7-i 
Jl=2U&*n3U 
J 2 = L "■ I S * M 3 'U' ! 3 
FQ-R K-=J1 TO J2 
Z(r .-^Z- K-I '-^Z4 
•y K / = Y - 1 K - 1 ; 

NE'-T K 

J3-;->I6*M3-K13+i 
2( J3;=Z(J5-l > 

•(UU ; = Y< J3-1 >-^Y4 
J4=J3rl 

J5=2*IS*M3+2*M3 
FOR K=J4 TO JS 
2 -IK ;- = 2^K'-l >-24 
Y(K > = Y(K-1 ) 

NEXT K 

IF n=N3 THEN 1010 

JS=J5+1 

Z( Jb ) = Z( JO-1 

r< J6 >=y< Jd- 1 !-n'4 

NEXT 17 

REM 



S3 



1 



1030 ASSIGN C^PatnZ TO "PRSl’" ' 

1031 INPUT "DG you WANT TO SKIP PREulNlriARr ’ |j 

1032 IF N«»“Y“ THEN 2290 
1340 REN 

1050 REN THIS SECTION COMPUTES Tn£ ZE-G OPFSET CORPECTICN 
iGbO REM FOR THE TR-^NSDUCERS . 

1070 REM ■ 

10S0 FRINTER is ! 

1090 PiTIhI “fOMPUTE THE ZERO OFFSET CORRECTION" 

1100 T-RiNT "Mi Oiccor-JMECT ALL TUBING TO THE TRANSDUCERS m i 
1110 PRINT " HIT COriTINlJE WHEN READv" 

1123 PAUSE 

1133 FOR 1=230 TO 204 

1143 G2=3 

1153 FOR J-=i. TO 50 

1 1 S0 OUTpijT 709 •."AIM," IT 1 " 

1170 ENTER 709,.kZ 
11 £3 G2=GZ-^K2 

1190 NEXT J 
1200 G5=G2/ 58 

1210 64( r -199 )='33 i 

1223 NE-T I 

1230 REM END OF LOOP FOR ZERO CORRECTION OF TRANSDUCERS 

1240 pen 

1250 REM THIS LOOP CALIBRATES THE TRANSDUCER AGAINST THE MANOMETER 
1260 REM 

1273 FOR 1=230 I'O 2?i 
1280 11=1-199 

!1'30 »^RINT •• TRr^N50UCER CALIBRATION, ENTER THE MANOHETFR PRES3URR I^i i'nPh-S O'" 
HATER" 

1300 PRINT 
1310 PRINT 

1320 PRINT "ENTER THE PRESSURE FOR TRANSDUCER NR . " , I i 

1330 INPUT H2 

1340 G0=0 

1353 FOR J=i TO 53 

1360 OUTPUT 739: "nl" - I ; "UTr' 

1570 ENTER 7QS;X 
1380 G0=G0+X 

1380 NEXT J 
1400 El=G0/50 

1410 A( r-199)=H2/<El-S4(I-iS9) ) 

1420 PRINTER IS 701 

1430 PRINT "A" ,A( I-1S9 > 

1431 PRINTER IS 1 

1432 PRINT "AM," = ",A(I-1S9) 

1440 INPUT "IS AI 1-199) SATISFACTORY rv/NCPNS 

1450 IF M'5 = "N" THEN 1293 

1460 PRINT 

1470 print 

:490 Ni:-'T I 
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iSC0 

1519 

1523 

1530 

1 54-9 

1550 

15B0 

1579 

1530 

1590 

iS00 

1GI0 

1520 

1G30 

1540 

IS 53 

iSS0 

IS70 

1533 

15 99 

1700 

1710 

1720 

1730 

17.:0 

1 ?5G 

1753 

1 770 

1730 

1790 

1S30 

ISIG 

1829 

1S30 

1840 

1850 

1360 

1870 

1300 

18S0 

1300 

1910 

1920 

1930 

1940 

1950 

i 980 

1970 

1930 

1990 

2000 



REn END OP LOOP FOR CnLIBRPTrMG TRhmSOL'CERS 
REM 

REM ENTER THE OAluES OF bTATfC AND TOTAL -'RE53URE OF THE FREESTPEAM 
PRINT •'INPUT '5TATIC FRES'T'URE OF FREETTREAM IN.HZOr' 

PP.INT -‘USE TRANSDUCER NR.! AND CHANr^hL f.jR.28'J- 
PRINT 

PRINT "HIT CGNTIKUE 1;h£N P •STATIC I? PRL'AEF'uV SET ;jp" 

PAUSE 
B S “ 0 

FOR J = 1 TO 5-0 

CUTP'JT 739 ; fil - ; 2S0 ; -UT 1 " 

ENTER 709 iXB 
BS =b6+ k 6 
NE'lT J 

PEN EMC LOOP FOR ACGUIRiNG STATIC PRESS. OF r REE2TREr,n 

P 9=56/ 53 

print 

PRINT 

REM 

PRINT '•INPUT TOTAL PRESSURE OF T.HE FREESTREAN, < [m.h 20 
PRINT “USE TR.ANSOUCEP NP.l AND CHA.N.N.EL NR'.200‘' 

PRINT ■‘•AIT CONTINUE UHEN P TOTAL PROPERff 'SETUP" 

PAUSE 

B7=0 

FOR J=i TO 50 

OUTPUT 70'3 ; “ A I " ; 200; “ OT 1 " 

ENTER 70S;X7 
E7=B?-r’T7 
n£>:T J 

REM END OF LOOP FOP fiCOUIRIfiG TOTAL PRESSURE OF FREESTREAM 

P7=B?/S0 

REM 

REM ENTER AMBIENT CONDI T ION'S 
INPUT “ENTER FAMB (IN. OF H;3i",Al 
REM 

REM LOOP FOR ACQUIRING TE-.MP OF FREESTREAM UIA rHERMOCOUFLE 
T=0 

FOR J=l TO 50 

OUTPUT 70S; ”AI“ ; 109; "UTl " 

ENTER 703: T1 
T=T+T1 
NEXT J 

REM END OF LOOP 
T2=T/50 

REM CONVERT VOLTAGE TO TEMFEPATUPE 
Ei=T2*1000 

T3=26.573*El-i.93637S'El '.:+ . 387a5*E T'3- . .25 1277 ‘E 1 ' 4 
REM 

REM CONVERSION TO SI UNITS 
T4=T3+275. 15 
A 1 'A 1 ‘3385.82 
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2S10 Rj=Ai/';2S?*T4 ,) 

2020 C = <P7-P9)<^2^S.7 
2030 REH FREESTREAM VElOCITY 
2040 U3 = ^ 2*C/R1 ,S 

2053 REi'1 

2060 REM BEFORE EACH OATA RUN ENSURE THAT P2=P3 i.E. THE PROBE IE BALANCED 
2070 PRINT "BALANCE F2 AND P3. " 

2080 PRINT "HIT CONTINUE WHEN ALL rUBIf'iS HAS BEEN CONNECTED" 

2090 PAUSE 

2100 REM LOOP FOR SAlANCIMB PROBE 

2110 FOR 1-201 TO 202 

2120 Gl=-0 

2130 FOR J=1 TO 53 

2U0 OUTPUT 70S; "AI" ; T ; "UTl" 

2150 ENTER 7G3;X 
2 160 G 1 -b 1 + 

2170 NEXT J 
2180 X-Gl/50 

2130 E'. 1-199 )=>; 

2200 Pi 1-139 '=A'- 1-139 )*E< 1-193 )-G4( I-I9S ) 

2210 NETT I 
2220 PRINT 
2230 PRINT 

2240 PRINT "P2=",P(2) 

2253 PRINT -P3-“,P<3) 

22S3 INPUT ‘'DOES P2-P3 (Y'-N/",N4 
2270 IF N$="N" THEN 2100 
2280 REri 

2233 REM ENTER THE LOOP FOR ACQUIRiNS EACH FRESSuF.E COMPoriNC CJEPF'I CENTS 
2330 REN AND COMPUTING TUO PRESSURE COEFFICIENTS AND TWO PRESSURES 

2301 PRINTER IS 701 

2302 PRIriT "Y 7 Mi M2 PI PAUG*- 

2303 PRINTER IS 1 
2310 REM 

2320 K3=0 

2321 !;2=M3*N3 

2330 FOR K=1 TO K2 

2331 K9=K 
2340 REM 
2390 WAIT 10 

2400 REM LOOP FOR ACQUIRING EACH PRESSURE 
2410 FOR 1-200 TO 204 
2420 Gl*0 
2430 REM 

2440 FOR J=1 TO 10 'ENTER THE DAS AND SAMPLE EACH PRESSURE IQ TIMES 
2450 OUTPUT 709 : “ AI" ; I ; “UT i " 

2460 ENTER 709; X 
2470 Gl-oi+X 
2480 NEXT J 
2490 REM 

2500 X-Gl/10 lAUERAGE THE VALUES 
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CD CD CD uli c:* CO CD 0 -‘ CO cu CO ro I.U CO -j vj - j nj ..I -.j -j -j o'j m cn rri r/:i cj m a* u'\ m ai u’l lTi li'i u t in ui i.n m 



Ei I-i9S * = X 

P ( I -* 1 99 ’=f^< I “ i 35 / '*‘E \ i - 1 'j‘J ;-b4 I - 1 99 * CUPPECT I OrJ FOP ZERO AND TRAnS. 

IF I-ZO0 THEN 

P1(K )-P( I~1S3 > 

END IF 
NEaT I 

PEP END OF LOOP FOR ACDUiRlNG RRE9SURE9 
^‘EN 

KEN CuNPUTE THE AvERnSE OF P:,-3,P^,Pb 



i - D 



FOR 
!JE=f" ( I .*'^'05 
NEXT i 

b o ( ^ b / ^ 

P£M END LOOP 
REN 

RE'1 CDriPL' i b : Ht UObFF I C i cMTS O'* VAN nND r J. TCM 



I LOCA-L ypw‘ CCbFFiL'ENT 
! LOCAL PITCH COEFFICEmT 



2350 
2 330 



L>l-P< 2 )-F( Z ■ 

02’=FR 4 )-P{ 5 ) 

DZ~P ( 1 ? —r 6 < K ) 

Ni(K )=Dl/03 
riZ' K '=02/03 
PFI.NTEP is 701 

print using 2750 : Y( K ; , 2; K ) , PU i k \ N2( K ^ Pi C K ) , PS( K ) 

inAGE NnD.C0.:x,nD[;.DD,2X ,riDD.0DDD,2Y , MOD , DDC D , 2 >! , NDQ . ODDD , 2X , MOD . DGDD 

R£M 

riEv.T K 

PEN 

FOR 1 = 1 TO KS 

OUTPUT 5p3th2;Y( I ) , 2 • I ;,M1- I ).M2( i ),Pli I ;,PS( I ; 

NEXT I 

nSSIGN OFathi TO » 

ASSIGN OPathC TO ♦ 

PRINT <-**.*-^r4444^444'* 

-•RINT 

PRIMT 'VELOCITY COfiPUTftTIOr-; VELOCITY COfiPONENTS' 

PRINT 
FP IilT 

PRINT “CATE OF RUN IS" JiS 
PRINT “TIME OF RUN I5“,N9 
PRINT 

PRINT “DErJSITY(KG/M-3 )•■ .P.l 

PRINT 'FREESTREAM VELOCITY (M/S;",U3 

PRINT “PAMBIENT(N/M- C ,A1 

PRINT ''TAMBIENT<C)“,T3 

ENO 
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ic' 

20 

3S 

40 

50 

S0 

■20 

30 

50 

100 

no 

120 

150 

140 

150 

ibO 

170 

ISO 

1 30 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

3G0 

310 

320 

330 

340 

350 

3S0 

370 

330 

390 

400 

4 10 

420 

450 

440 

450 

460 

470 

480 

4 J0 

500 



REn PRCGRAM vEL 

REM THIS PROGRAM COMPUTES THE COMPONEMTS OF UElOCITY FROM 

REM THE RAW DATA COLLECTED UIA THE PROGRAM PRSACQ 

REM 

REM AUTHOR OAL'E EUAME , .lANUARV . FEBRUARY , 13-27 
REM 

REM VARIABLE NrtMtS 
REM 

REM FS=AyERAOE PRESSURE 
REM 

REM rn -MEASURED YAD CGEFFICIEMT 
REM M2=M£A3URED PITCH COEFFICIENT 
REM m 3=MEA5UPE0 TOTAL COEFFICIENT 
REM M4=n£ASUR£D STATIC COEFFICIENT 
REM 

HEM l'=THE TCmAL 'v'ELOCITY AT THE P-OSE TIP 

REM Ui'THE X COMPONENT OF UELOCITV 

REM U2»ThE Y COMPONENT OF VELOCITY 

REM UZ=THE 2 COMPONENT OF DELOCITV 

REM 

REM Xl-THE PITCH ANGLE FROM THE CALIBRATION DATA 
REM Y-THE YAW ANGLE FROM THE CALIBRATION DATA 
REM 

REM X4-THE PITCH ANGLE OF THE FLOW 
REM 

REM y=THE VERTICAL LOCATION OF THE PROBE TIP 
REM Z=THE SFANWISE LOCATION DF THE PROBE TIP 
REM Y3=TH£ VAN ANGLE OF THE FLOW 
REM 

DIM C 1 < 77 ) , C2 ( 77 ) , C3< 77 i , C4 ! 77 ) , MU 80C > . M2 ; 200 > 

D 1 M A 1 ( 77 ) . Y V 77 ) . V 1 ( 80S ) , P I ( SO0 ) , Pd >: SCO . P S ■; 200 ! 

DIM 20< 300 ) ,0<; 600 ) ,U1 i S00 , L'2( 800 ) , 031 830 • 

DIM Z ( 800 ) . FH S03 ) , 04 ( 800 ) , 05 ( S30 ) , UG t 800 , '■/? ( 300 ) 

R E M 

INPUT "WHAT IS THE DENSITY FOR THIS RUN'. KG/fr3 i " . R 1 
REM 

REM READ CALIBRATION DATA INTO COMPUTER MEMORY 
REM 

ASSIGN QPathl TO "CAL" 

F'OR 1 = 1 TO 77 

ENTER §Pathl;Xl(I ),Y> I ^C1( I >.C2k I .v,C3':i 1,C4< I ; 

NEXT I 

REM READ RAW DATA INTO COMPUTER riEMORt 
ASSIGN @Path2 TO "PR52'' 

FOR 1 = 1 TO 8'00 

ENTER ■@Path2:Yl(I ),2H ),.M1( I >.M2( I > , P 1 d / , PB< I '■ 

■NEXT I 
REM 

CREATE GOAT '•U£L2".2 25 
ASSIGN OPathZ TO "UEL2' 
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51£ FOR K = 1 '0 <?08 
520 K9=K 

530 >'1;k )=i 1(K )-3.5-i 

540 Z;K ) = -2'K )‘:;.£4 



550 
aSO 
S^O 
530 
590 
500 
□ 10 
s:o 

630 

540 

550 

6SQ 

570 

SS0 

590 

700 

713 

7:0 

750 

743 

750 

750 

77G 

7S0 

7 30 

303 

310 

530 

G3G 

840 

950 

3S0 

•373 

930 

893 

330 

913 

320 

330 

S40 

353 

3B0 

970 

993 

330 

IG03 

i 0 1 3 



REn 

REr-l COnFUTE THE FIRST APPRO;-; In, -VT 10m OE ariOLE 
B3--. 1538S324 

B4=— .3555£718 ' COEFFICEMTS ANO POLVr-JOrilAnRL ORE 

0-5=. 13550589 ! rROM FORTRnM PROGRATl "POLV AMO 

Gb-.3C401553S I ARE OCMPUTED FRCM THE CP vfto OS. 

C7=. 033000069 ! Y.AO DATA FROTi THE CALIBRATION 

C9=-.030OSB67 

■:'3=B3^S4^Miv:-H )+C5>rilCK )"2 ^CS♦^^1U^ :'"3+C7-i1l(K !-4+C9 ‘Ml ( K r 5 
»EM 

FEM GO TO the SUB-ROUTINE FOR PITCH ANGLE InTERPOLATIOM 

C-j5U0 14 00 

HEM 

RE" GO TO THE SUB-ROUTINE FOR VAy ANGLE INTERFOLAT i ON 

dU 3UH i U j3 

R£r, 

HEM GO TO THE SUB-ROUTINE FOR CP TOTAL AND CP STATIC INTERPOLATION 

rObijP 135.11 



p ;□ M 

REM THIS SECTION CO.MrUTES THE MAGNITUDE OF TH£ LOCAL UELOCITY 

PEN vEGTG.R AT THE PROBE TIP 

PEN 

PKK :-=Pl<. K !'Z4 8.3 ICONMERT PI TO N/N'E 
P6vi- .!-PS(K ;-249.3 iCONUERT PCS) TO N/M-2 
.HEN calculate •STATIC AND TOTAL PRESSURE 
r= = P5;.C '--NP* PKK )-PG(K :■ > 

FS = -P8 lOO.NE BECAUSE INTERPOLATION RESULTS IM NEGATIUE UALUE FOP PSTAl 

? ; K .'=F 1 ■■ K • -N3 c P 1 < N )-P6 T k ? '< 

U4! r. 2- Ri P(K >-P3 ) )•■ .5 



REN 

REM COMPUTE THE COMPONENTS OF UELOCIT'i 
CE6 

USfK ;.-=U4(K )*C0S1 X4 )*C0SCr3 ) 

■US'C: )=U4(K )*SIN1X4 ) 

■U7! K >=U4{ K )*CQS< X4 ) *3IN< Y3 ) 

NEXT K 

REM END OF DATA REDUCTION LOOP 
REM RE-ORDER DATA 
11 = 1 
12 = 40 

FOR 1=1 TO 721 STEF 80 
Jl = Ii 
J2 = I2 

FOR K=1 TO 40 
7.0! Ji ) = :•' J2 ) 

P3( Jl )=P( J2 : 

Jl )-U4v jc •■ 
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1020 Ui Ji ='v'5ij;: i 
1030 y2< J] ,'-06( J2 ) 

1040 y3( J1 J2 ) 

1050 Ji=Jl+i 
1068 J2=J2-i 

1070 MEXT K 
1880 11=11+80 
1090 12=12+38 

1100 rJEXT I 
1110 13=41 

1128 FOR I =4 1 TQ 808 SIEF- 08 
1138 K=I3 
1140 FOR J=1 TO 40 
1150 Z8(K) = Z>'K) 

11G8 P91K /=F-:k ) 

1 J 7D \/{ ' = '•■' 4 i R > 

list) OliK) = OS':K} 

1190 •.■'3(K :'=--U6(K ) 

1208 F'3(K '=y7{K ) 

1218 K=K+1 

1228 NEXT J 
1230 13=13+30 

1248 NEXT T 
1250 FOR K = 1 TO ,-;2 

1 268 OLi TPUT fpPa t h3 ; 71 1 !< ) , Z0-: K '■ . F3' F ' , I K . 
1270 NEXT K 

12S8 ASSIGN @Pat hi TO ‘ 

1230 ASSIGN @Path2 TO » 

1300 ASSIGN @Path3 TO » 

1510 PRINTER IS 1 

1320 PRINT “Y(Cr1: Z'CM; PTuT (fj.rr 

vzai/S)" 

i.i38 FOR 1 = 1 TO K 9 

1340 PRINT USING 1350 ; Y i (I > . Z0( I ; , FS.; I > . 

1350 Ii-InGE MDD.DD , 2X ,nOO. DO , 2X .NOOn OOP 
.ODD 

1350 NEXT I 

1370 PEN 

13S0 STOP 

1390 REM 



I 



vl‘K ).02; K K > 

3' '-■•'N/S) OX(M/S) 0X1 M/; 

1 ! , VI ; I ' . 02 ■: I > , vj; i ) 

. MOD . DUD , 2X , non . ODD , 2x , MOD . ODD . 2 



1400 REM PITCH ANGLE INTERPQi riT ION 
1410 REM 

1420 FOR 1=1 TO 7? 

1430 IF-Y3<Y(I) THEM 1450 * SEARCH THRU CALieRATION DATA TO FIND 

^‘i‘^L+S WHICH BRACKET THE APPROX YAW 
1450 K1=I lAMGiF 

1460 K2=I-7 

1470 K5=I-i 



1480 K4=I+S 



1430 FOP j=K2 TO K3 

1500 IF M2 : K .1 C2>' J > THEM .(528 



'FIND THE PITCH ANGLE WHICH CORPESPONOS TO 
'The me {■'< s u e o >■- 1 t c h c; o e f f . 
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lEiO 
1520 
iS30 
1540 
15SG 
iSSG 
1S7S 
iS8Q 
1530 
!S53 
IclC 
3 620 
1330 

1O40 

IS 50 
1330 
1370 
1 560 
IBSO 
1700 
1710 
1720 
1730 
1743 
1750 
175 3 
1770 
176? 
i':’93 
1803 
1310 
1323 
is 33 
1840 
1350 
1853 
1370 
iSSO 
1390 
1908 
1913 
1920 
1933 
1940 
1950 
1953 
1970 
19S0 
i'3S0 
2300 



NEXT J 

Z1 = >C2U 1- i)/>C2(J !--C2(J -i -) 

X2»Kli'J)-2i»<XHJ>-Vi<J-i )> ; BRACK STING FITCH nfia E 

FOP J---r i TO K4 

IF MZa;? C2<J) THEN 1573 

NEXT J 



22= C2'. J '-02', K :■ ).-N C2l J ^-C2'. J- 1 :■ 1 

:'.Z=XM J .'-Z2‘iXl( J '.-XK J-I ) : IBRACKETING PITCH 0N5LE 

Z3 = ^' VI S4 ,’-V3 )/• ( V< Ki ; - Y( K3 ) ) 

x4 = x3-Z3-< X3-X2 ) iPITCH ANGLE OF THE FLOW 

RETURN 

REM 

REN VAN angle interpolation 

REM 

FOR r=l 10 77 
IF VZ’i -. (!'■ THEN IB80 
NEXT I 
i-'2 = I-? 

H3=I-i 



POP. J=HZ 10 H3 

IF X4 /.i;J) THEN 1730 

rJEXT J 

F1 = J 

F2=J-l 

FZ-j+£ 



IXIMO’ YAW ANFLiiS 3HiCH 
iCGRRESFONu TO TaE MEAS-UPEij 
iCP YAN ON PITCH ANGLES OHICH 
'BRACKET X4 



F4 = J-^7 



Z4-(Ci ' 7 1 ; -Ml< K ) )/< CKFi )-Ci<F4 ? ) 
v 4 = V‘;F1 )-Z-;‘>(Y(Fl l-Y(F4)> 



: 5 = ( C 1 • F Z ) - M 1( K ) > / < C 1 1 F 2 ) - C H 3 ) ; 
Y5 = v'^F2 )-25*( 7<F2 )-YN F3 ) ) 
rb=^>:i<Fi )-x4 ) ■'! XKFi )->:i;F2)i 



Y3=Y4-Zfa'i V4-Y5 i 

RETURN 

HEM 

REM INTERPOLATION FOR LOCAL COEFFFICIENTS OF STATIC AND TOTAL PRESSURE 

23=0 

Z2=0 

FOR I=i TO 77 

IF Y3-Y(I> THEN ISIO 

NEXT I 

Fi = I 

F2=I-7 

F3=I-1 

F4=I+5 

FOR J=F2 TO F5 

IF X4 .X1(J) THEN 1980 

NEXT J 

F5=J 

Fb=J-i 

F7 = Jr'7 
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rj rj 



Z010 F2-J+6 

20Z0 Z3=( X U Fb ;■ -X4 / ■' ( >; 1 ■; F6 .'-X l < FS ) ,) 

2030 N4--C4< FG )-23* ( C4! Fb )'-C4C Fb ) 

2340 N3=C3: F6 '-Z3' ( C3< F6 > -C3; FS :• > 

053 L4=C-(F3 )-23'<C4(F3 )-C4< F7 ) ) 

0S0 L3=C3( FB ;-Z3^ ^ C3; F8 ‘"-C3-: F7 ) ) 

2073 Z2-( Y(FS ■-> 3 Y< FG )-V;F8 ) ) 

2080 !-14-N4-22-' N4 -L4 ,» ! LOCAL STATIC COEPFICENT 

20S0 M3-=N3-Z2»(N3-L3 ) i LOCAL TOTAL COEFFICImT 

2100 RETORN 
2110 STOP 
2120 SriD 
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I 

4 

20 

30 

50 

50 

70 

80 

50 

100 

i 10 

!2G 

130 

140 

150 

150 

17G 

180 

13G 

200 

210 

220 

221 

223 

22^ 

225 

22 £ 

230 

Z'iS 

250 

260 

273 

2S0 

2^0 

300 

301 

302 

310 

323 

321 

330 

340 

341 

350 

3S0 

370 



REM PRGGRRM PLOT 

REM THIS PROGRAM MAKES CONTOUR PLOTS FOR STREAHVI3E 

REM UELOCITY AANO CAN BE MODIFIED TO MAKE COHrO'JFS PLOTS - 

REM FOR TOTAL PRESSURE 

DIM V i (• 800 J , Z01 800 ),P9( 800 > . v'. 300 ' . Ui 1 300 ' . v'2 ' 808 1 . L'3< 800 ■ 
DIM A<S00) 

ASSIGN •SPethi TO “UELCr' 

POR I---! TO 800 

ENTER OFa thl ; 'i 1 v I J . 20; I :■ , p‘3( I 5 . vi I T, 0; ' I ; , v2 • I , ''.■'O [ • 

Nt’T ! I 
'.•3=V\ 1 > 
i ;■ 

FOR 1=2 TO 800 

IF U1-. I :'>U9 THEN 03=01': I 

IF y-vI)U3 THEN U3=yi<I) 

NEaT 1 

D = : U3-US )/I0 
FOR 1=1 TO 500 
Z3='JS<-D 

IF Vi ■: I ; ;Z3 THEN Xi 1 ,'=G 
FOR J=i TO 9 
2 i =v : J + 1 ! "D 
Z2='.'Sp i 

IF UliD'Zl AND Uiai'22 THEN :-:U ) = J 
^'L^’T J 
ME>:T I 

FOR I-=l TO 300 

IF X'll KS THEN GOTO 226 

ri“Ua=i 10-0 ) 

Z2=US-i-: 9.7*0 1 

IF uia-v:i and uiuiTZZ then >;(I)=0 

NEXT I 

GINIT 

'PLOTTER IS 735. ’‘HF-GL" 

PLOTTER IS CRT, “INTERNAL" 

8F'APHIC:3 OM 
CSIZE 4.5, .65 
MOUE 39,10 

LABEL “EMBEDDED UORTEX 28 M/S " 

MQL'E 43.5,7 

LABEL "75% FILM COOLING" 

MODE 41,85 

LABEL “STREAMWISE UElQCITY" 

MODE SB, 17 
CSI2E 3.0 
LABEL "Z CM " 

MQUE 5.5S 
CSIZE 3.0 
Label t=“T CM 
FOP :=-i TG G 

LA6EL LateUlI.n 
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I33i3 


li£XT I 




390 


VIEWPORT 15, 114, 25. 79 




400 


SHOW -iG,S,0,l-'l 




410 


FROi-lE 




420 


AXES .5, .5. 0.0. 2, 2 




430 


AaE3 .5. .5,-16,12.2,2 




4 40 


AXES .5. .5.6. 12,2,2 




4 41 


C5IZE 2.3..B5 




44 2 


r'OUE -1 4. 7, 10.55 




“ --r ^ 


LASEl. ■■VELOCITv RrAilG£5" 




4 4 b 


A=10. 1 




■i-iS 


FOR K=3 rO 3 




4 47 


S = A-K » . 35 




448 


lioue -ii.T.e 




4 48 


22---VSfK -0 




458 


21-V-if' t>i wo 




451 


IF 5=9 then ZJ=V84k'D+.?*D 




452 


label UEILG ‘■D,2;<,OD.O,X,DO.D" 


.K,Z2,2i 


453 


riExT K 




454 


nOVE -14.3,B-.35 




455 


LABEL US I rlG " D . 2X , 00 . D , X , DD . D “ 


; 0 . 2 1 , V9 


45c 


C5IZE 1.3,. 72 




457 


FOP 1=^1 TO SOO 




■159 


r-lOl'E JChI/.YI'I) 




453 


IF Yi-a L>5.0 ,AriD 201 IK-i0. THEll GOTO 451 


463 


LABEL USING “0“:X<i / 




45 i 


NEXT I 




4S2 


CLIP OFF 




455 


CSIZE 2.2. ,55 




4? 8 


FOR I = TO 11.8 STEF 2 




/■"'W 


HOVL -17.5,1 




430 


Il-If .2 




53S 


LABEL US If, '6 •‘^.DD.CW ; 11 




510 


NEXT I 




5Z3 


FOR J=-1?.0 TO 5.0 STEP 2 




SZ?j 


nOVE J,-.7 




540 


J1-J4 1 .0 




550 


LABEL US I MG " t , MOD . C “ : J 1 




560 


NEXT J 




S70 


END 
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12 

20 

30 

40 

50 

BO 

70 

30 

30 

100 

110 

120 

130 

U0 

150 

ISO 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

230 

233 

330 

310 

323 

330 

343 

353 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

4E0 

470 

480 

4S0 

500 

c; ! 



REn PROGRAM UECTOR 

REM THIS PROGRAM USE DATA -ROM ULLC DAI A FIlES 
REM TO PLOT SECOf-iOnRY RLOU v'EOTOPS 

DIM VI ; 830 ' . Z0( 800 ' ,P9; 800 • ,U! 902 ), v'l ( 900 J ,'’ 2 < S03 ■>.'!}( 8G0 '• 
DIM X(303),V(5),Z;5'- 
ASSIGN L3Pathl TO "vELCin" 

FOR 1=1 TO S0O 

EfiTER gPathI ; VI v I ). Z0( I ) ,P3' I ) ,U i • ,:,M •: i ). v2- I I ' 

MEX f 1 
GINIT 

Plotter is 705,-hpsl- 

I PLOTTER IS CRT , " INTERNAL " 

GRAPHICS ON 
CSI7E 4.5..G5 
MOVE 39, IG 

LABEL "EMBEDDED U0K'E:< l1 N 3 " 

MOVE 43.5,7 

label "OS?: FILM OOOLUiG" 

MOVE 35.4. 85 

LABEL “SECONDARY FLC.W VECrOfiS'' 

MOVE 52.1? 

CSI2E 3.0 
LABEL "Z CM “ 

MOVE 5. 65 
C5I2E 5.0 
LaLel5="V CM " 

-OR 1=1 TO 5 

L.A3EL Labe 15 [ 1 , I ! 

NEXT I 

VIEWPORT 15,114.25,79 
SliOw ~ 1 6 , 5 , 0 , i 2 
FRAME 

AXES .5, .5. 3, 0.2, 2 
Axes .5, .b,~l6, 12. 2, 2 
AXES .5, .5,6,12.2,2 
CSIZE 1.3, .72 
Sl=.21 
A 1=0. 

FOR K=7Si TO 800 
Al=AlxVUK ) 

NEXT K 
A1=A1*SI.''40. 

MOVE -14.5,10.0 
Y< 1 )=10.S 
Y<2 )=10.B 
2< 1 >=-14.5 
Z^2)=-14.5+nl‘.l 
FOR J=1 TO 2 
Plot z<J),v<j; 

NEXT J 
PFNVi'" 
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520 

530 

540 

550 

5G0 

570 

550 

5SS 

503 

G10 

G20 

■530 

5 40 

550 

B60 

S70 

580 

593 

700 

7 13 

720 

738 

740 

750 

750 

770 

788 

700 

S00 

810 

■820 

S30 

340 

35S 

858 

a 7-3 

8S0 

890 

933 

910 

920 

330 

940 

950 

9B0 

370 

988 

S30 

1000 

1313 



M0y£ - 1 5 . ,3.5 

CSIZE 2.0, .55 

LrieEL 'TJ.l FS 'JEL0GIT7" 

FOR 1=1 TO 303 
Ys 1 ) = V1': I ) 

2( 1 !-20; I ) 
r ( 2 ) = V1 i I i-»SI‘U2\ I ) 

4'. 2 ! •-'•Sl*03( I ' 

I 

! BOS I MESS TO MAKE ARROWS 
! 

Di-' V‘ 2 )-Y( i ) ) 

D2 = ( Z{2 )-I< i ) 

DEG 

A3=--ATN{D2/01 ) 

A 4 = ■i 5 ” A 3 
£3 = 3, 

rr CIO AND D2;0 then E3=1.0 

IF £3=1.0 AND A37-45. THEN A4=.A3+i 35 

D4=i ( 2; 2 )-2( 1 ) r'2 + ': y< 2 ,>-Y( i ) r'Z r" .5 

L’4=.2G‘04 

06= i .3 

03=1.0 

IF A-; 4 5 THEN D3=-l.0 
E1=S. 

E2=8. 

IF Dl:-0 AND O2<0 THEN 09=-l.e 
IF Di:0 AND D2--.3 THEN DS=-1.0 
IF 01 0 AND D2<3 THEN D8=-i.0 
IF Gi-OO AND D2<0 THEN D9=-I.-3 
IF Di'-73 AND 02,- 0 THEM D8 = -1.0 
IF O1'0 AND A3>45 THEN 03=1.0 
IF A3;--i5. AND A5<0. THEN D8=1.8 
IF D1>0 AND D2<8 THEN El =1.0 
IF £1=1.3 AND A3:-45. THEN 03=1.0 
IF OK0 AND D2<0 THEN E2=1.0 
IF E2 = 1.0 AND A3 .45. THEM D3=1.3 
V(3) = Y(2 )-D44C05<n4)*D8 
Zi3)=Z<2)-D4*9IN(A4)-D9 
Y( 4 ) = Y(2 > 

Z<4) = Z(2) 

A5=A4-90 

IF Ol>0 AND D2<0 THEN D8=1.0 
IF D1<0 AND D2<0 THEN D3=1.0 
Y(S) = Y(2>+D4*CQS(AS)‘09 
Z(5 )=Z(2 )+D4«SINl A5 )^D8 

I 

I END OF BUSINESS TO MAKE ARROWS 

I 

MODE I'Ol i I } 
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1030 


FOR J-l TO G 




1040 


PLOT Z( j J • 




1853 


NEXT J 




IC'BO 


PcNUP 




1 870 


NEXT I 




1030 


Clip off 




1 0 9 vj 


CSIZE 2 ,0, .95 




i 1^0 


PQP TO i 


1 . £ £TEP 2 


1110 


riO'JE 




i : C'iy 


J 1-I+.2 




I 1 .jG 


' R E L U S ^ N G 


' It , 00 . D " 1 ! 1 


i 14G 


NEXT I 




il50 


FGR J — 17.0 TO 


b.0 STEF 2 


IISG 


NOOE J,-.7 




1170 


J ]. ■= J i- 1 ^ 0 




1186 


LABEL USING 


“ ? , i'lDD . 0 . J i 


1190 


NEXT J 




1200 


END 
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APPENDIX C 

UNCERTAINTY ANALYSIS 



Uncertainty analysis was performed using the method originally attributed to 
Kline and McCormick,. [Ref. 13]. Let 6 j. be the uncertainty in the result and 82 . 
. ., be the uncertainties associated with each independent variable. The uncertainty 
in the result can be expressed as 

\ = [I((5R/5x.) X (6.))]1/2 (eqn C.l) 

To determine the uncertainty of the pitch angle, a, a straight line approximation 
was made for a 

a = m + b(Cppi^(,j^) (eqn C.2) 

The following independent variables were determined: = ±1.97 <>, = ± 

2.32®, and = ± 0.04637. The uncertainty of a, was determined to be ± 2.36®. 
The high uncertainy of the pitch angle is caused by the probe being highly sensitive in 
the pitch plane. 

The uncertainty of the yaw angle, P, is determined in a similar manner. Once 
again, a straight line approximation is made. 

P = m + b(Cpyg^) (eqn C.3) 

In this case, the uncertainties of the independent variables were found to be 6 ^^ 
= ± 0.954, = 0.5924, and 6 ^-p = ± 0.058. From these values the 

uncertainty of P was calculated to be ± f.29 ® 
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